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GENERALIZED RING STIFFNESS MATRIX FOR 


RING-STIFFENED SHELLS OF REVOLUTION 

By George E. Weeks and Joseph E. Walz 
Langley Research Center 

SUMMARY 

A generalized set of nonlinear boundary conditions and ring stiffness coefficients 
are derived for a shell of revolution with an elastic ring of arbitrary cross section 
attached to its boundary in an arbitrary manner. The effects of shear deformation, inter- 
nal pressure, eccentricity, restraint of warping, torsion, out-of-plane bending, and pre- 
stress of the ring are included in the analysis. The boundary conditions are applicable to 
shells with an axisymmetric initial equilibrium state of stress and deformation which are 
subject to small asymmetric deflections away from the initial state due to either bifurca- 
tion buckling or the addition of an asymmetric load. The results are presented in a form 
which can be easily introduced into a computer program for shells of revolution. 

INTRODUCTION 

Shells of revolution with rings attached at the boundaries are common structural 
elements in aerospace vehicles. A typical example is the envelope or "aeroshell" of a 
spacecraft designed for entry into a planetary atmosphere. The analysis and design of 
such structural elements require a knowledge of the edge restraint or boundary conditions 
on the shell provided by an elastic ring of arbitrary cross section. However, because the 
stiffness characteristics of an arbitrary ring are very complex and lead to complicated 
boundary conditions for the shell to which the ring is attached, analytical procedures for 
taking into account ring behavior are usually based on approximate and many times inad- 
equate theory. 

At the time this study was initiated, the most comprehensive treatment of ring 
boundary conditions in the literature was given by Cohen who presented in reference 1 the 
boundary conditions for an arbitrary shell of revolution with an elastic ring of arbitrary 
cross section attached to its boundary. Effects of prestress, shear deformation, and 
restraint of warping in the ring were neglected. Later in reference 2, Cohen included 
some prestress effects, but still, his analysis is only applicable to a ring whose shear 
center and centroid coincide, and in addition, the normal to the shell at the point of 



attachment must pass through the ring centroid. The present theory does not have these 
restrictions and has been developed independently of the theory of reference 1. 

The purpose of this paper is (1) to present the assumptions and theoretical develop- 
ment of a ring theory for the nonlinear boundary conditions of a general shell of revolu- 
tion with an elastic ring of arbitrary cross section attached to its boundary, (2) to present 
the resulting equations (ring stiffness coefficients) in a form which can be easily incorpo- 
rated into a computer program for shells of revolution, and (3) to demonstrate the use of 
the equations with selected results which show the influence of ring flexibility and ring 
eccentricity. The ring theory is general and includes shear deformation, restraint of 
warping, out-of-plane bending, eccentricity, internal pressure, torsion, and axisymmetric 
prestress of the ring. (For stress analysis, ring prestress refers to any initial internal 
load in the ring with no external load on the shell. For stability analysis, ring prestress 
also includes internal loading in the ring prior to buckling due to the prebuckling deforma- 
tion of the shell from applied loads.) 

The ring theory is derived on the basis of the total potential energy of the ring. The 
energy of restraint of warping and torsion are included in an approximate manner similar 
to that of straight -beam theory whereas the energy of bending and extension are consid- 
ered rigorously. The results are also applicable, with certain noted exceptions, to the 
special but important case of a pressurized ring of closed cross section. 

SYMBOLS 

The units used for the physical quantities defined in this paper are given both in the 
U.S. Customary Units and in the International System of Units (SI) (ref. 3). Conversion 
factors pertinent to the present analysis are presented in appendix A. 

A,B,C,D,F functions of displacement variables (see eqs. (3)) 

A r cross-sectional area of ring 

A r area enclosed by ring of closed cross section 

aj geometrical expressions defined by equations (B3); i = 1,2, 3, 4, 5 

bj constants defined by equations (B5) and (C14) which relate ring variables to 

shell variables; i = 1,2,. . .,21 

E Young's modulus of ring and shell 
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circumferential strain of centroidal axis of ring 


e x ,e z thickness-compression strains at any point in ring cross section 

e X z> e yz> e yx shear strains at any point in ring cross section 

ey circumferential strain at any point in ring cross section 

F r ,Fg,F a middle-surface stress resultants at shell boundary in cylindrical 

coordinate system (see fig. 1) 

G shear modulus 

G^j elements of 4 x 4 stiffness matrix 

% ,A>^ij,B ring stiffness coefficients defined by equations (El) and (E2), respectively 


G ij " G ij,A 
h 

^X^Z^XZ 


J 

^i 


^x>M z 

M v 

m v 

N 


+ G ij,B 

thickness of shell wall 

moments and product of inertia of ring cross section with respect to ring 
coordinate system 

torsional constant of ring cross section 

elements of 4 x 1 column matrix in boundary-condition equation identified by 
external loadings 

bending moments in ring 

meridional bending-moment resultant in shell 

nondimensional meridional bending-moment resultant in shell, M^yA^h^ 

hoop force in ring 

middle-surface stress resultants in shell with respect to intrinsic coordinate 
system (see fig. FI) 
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n 


Fourier index 


P 

Pi 

Qi 

<*n 

a*. 

q 

R e 

r 

r o 

r s> r c 

S 

s 

s max 

T 

*V 

U f 

U P 

U r 


internal pressure in ring 

geometrical expressions defined by equations (D2); i = 1,2,. . .,23 
generalized forces defined by equations (Dl); i = 1,2, 3, 4, 5, 6 
Fourier coefficient in expansion for normal pressure on conical shell 
constant radial line loading 

radius of curvature of shell which generates shell circumference 

radius of middle surface of shell at any point along shell meridian (see fig. 1) 

radius of middle surface of conical shell at large end 

radius of line of shear centers and line of centroids of ring cross section, 
respectively 

GA r r§ 

shear-stiffness parameter, — 

2EI X 

meridional distance measured from large end of conical shell 
total meridional distance of conical shell 

EA. r^ 

extensional-stiffness parameter, — - 

EI X 

/ \ / 

nondimensional meridional stress resultant, / q n r Q 

potential energy of external forces applied to ring 
potential energy of pressurizing gas 
total strain energy of ring 
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u l> u 2> u 3 strain energy of ring bending and extension, restraint of warping, and 

shear, respectively 

u,v,w axial, circumferential, and radial displacements, respectively, of shell 

middle surface in cylindrical coordinate system (see fig. 1) 

axial, circumferential, and radial displacements, respectively, of any 
point in ring cross section in cylindrical coordinate system 

axial, circumferential, and radial displacements, respectively, of any 
point in shell wall in cylindrical coordinate system (see fig. B2) 

axial, circumferential, and radial displacements, respectively, of shear 
center of ring in cylindrical coordinate system (see fig. Bl) 

u,w meridional and normal displacements, respectively, of shell middle surface 

in intrinsic coordinate system (see fig. FI) 

w^ nondimensional normal displacement, w^ 

q n r § 

x,y,z ring coordinate system (see fig. Bl) 

x c ,z c axial and radial distance, respectively, from ring shear center to centroid 

(see fig. B3) 

x,z axial and radial distance, respectively, from ring shear center to point of 

attachment of ring to shell (see fig. B3) 

z o eccentricity of ring centroid normal to shell middle surface measured 

positively along outward normal 

z sh normal distance from shell middle surface in intrinsic coordinate system to 

point of attachment of ring to shell 

z distance of shell displacements u p , v p , and w p from shell middle surface 

a angle which meridian line of shell makes with radial direction (see fig. FI) 

/3 rotation of ring cross section about axis of ring shear center (see fig. Bl) 


u a> v a> w a 

u p> v P’ w P 

u s» v s»w s 
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r 


warping constant of straight -beam theory 


y ,y Timoshenko-type ring transverse shear strains 

yx yz 

6 operator which indicates variation of variable or function 

£ angle of twist of ring cross section with respect to its shear center 

77 meridional coordinate 

77 j quantities defined by equations (C3b) to (C3f); i = 1,2, 3, 4, 5 

6 circumferential coordinate 



£ = C + A' 

<p rotation of shell appearing in equation (B2) 

Subscripts: 

A initial equilibrium state 

B incremental equilibrium state 

Superscript: 

(n) Fourier coefficient of indicated variable 

A prime indicates a derivative with respect to the circumferential coordinate 6. 

ANALYSIS 

Statement of the Problem 

The most general boundary conditions for a shell of revolution supported by some 
generalized unloaded elastic medium at its boundaries can be written in the following con- 
venient matrix form (see, for instance, ref. 1 ): 
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The Gjj are stiffness coefficients which define the stiffness characteristics of the sup- 
porting medium; u, v, w, and (p are the middle-surface shell variables and F a , F$, 
and F r are stress resultants at the shell 
boundary in the axial, circumferential, and 
radial directions, respectively. The quan- 
tity is the moment per unit length 

imposed on the shell boundary by the elas- 
tic support which produces a rotation about 
a tangent to the circumference. Here, all 
terms are referred to the cylindrical coor- 
dinate system shown in figure 1 and the 
positive sign in equations (1) applies if the 
ring is connected to the positive (terminal) 
end of the shell. 

The column vector with elements ^ 

.. , .... , . , . , Figure 1.- Stress resultants, moment resultant, displacements, 

is nonzero if additional external middle- and rotation of she n a t boundary. 

surface forces on the boundary are specified. 

The effects of external loads on the ring itself are not considered in the present paper. 
However, the method for taking these effects into account is derived in detail in refer- 
ence 4. 

The basic problem is to determine the stiffness coefficients Gjj of an elastic ring 
of arbitrary cross section attached in an arbitrary manner to the boundary of an arbitrary 
shell of revolution so that the boundary conditions (eqs. (1)) for such a shell can be 
defined. The procedure used to determine the desired coefficients will be based on the 
variation of the total potential energy of the ring. 

Potential Energy of the Ring 

To formulate the potential energy of the ring, the displacements of any point in the 
cross section of the ring are referred to displacements of the shear center. The 
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displacements of the shell through the thickness are then determined and the ring dis- 
placements are expressed in terms of the shell middle-surface displacements in the 
cylindrical coordinate system (fig. 1 ) by requiring compatibility of displacements and 
rotations at the line of attachment of the ring to the shell and by assuming that shell dis- 
placements vary linearly through the thickness of the shell. These kinematic relation- 
ships are derived in appendix B. The final expressions for the displacements of a general 
point in the ring cross section u a , v a , and w a in terms of the shell middle-surface 
displacements u, v, and w; rotation 0 ; and ring transverse shear strains yy Z and 

r yx are 

u a (0,z) = B(0) + (z - z c )( p(6) 


Z - Z n X - x r 

v a ( 0 ,x,z) = C(0) - D(0) - F(0)S 


( 2 ) 


w a (e,x) = a(9) - (x - x^j(p(e) 


where 




A(0) = -w(0) + b i<p(6) 

B(6) = -u(0) - b2 0(0) 

C(0) = b3v(0) + b5w f (0) + biou f (0) + b 70 r ( 6 >) + bi5y yz (0) + bi7y yx (0) 

D(6) = -r s jb 4 v(6>) + b0W'(0) + t>i]u'(0) + b 8 4>’(0) + bieryz^ 6 *) + bl^yx^ 61 ) 
F(0) = -I , s [b 12 u , (0) + bg <p'(e) + bi8y yx (0) 




(3) 


J 


and where r s is the radius of the line of shear centers of the ring cross section and the 
coefficients bi are given b y equations (B5) in appendix B. 


The total potential energjr is composed of the internal strain energjr U r of the 
ring, the potential energy 7 Uf of the forces acting on the ring, and, for pressurized rings 
of closed cross section, the potential energjr Up of the pressurizing gas. 

Strain energy - The strain energjr for a circular ring, including the energjr of 
bending, extension, transverse shear, and the approximate energy of torsion and restraint 
of warping, is derived in appendix C in terms of the shell variables and the ring trans- 
verse shear strains. The result is 
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( 4 ) 


TT r 2 *[ Ne~ o «, e n 

Ur = fio ( 2 ^ 2 - Are o-^¥ Mj ‘ + -E s ’ , i + -Er’ , 2i rcdf ’ 


Er 

2r s r c J 0 


E 

^277 


E 


M z 

~e 


+ »2^") 2 r c de + ^|2 + r^r c do 


¥-^O b2o,# " +b2iu, ^ rcde 


where r c is the radius of the line of centroids of the ring cross section, the quantities 
e 0 and are given by equations (C3) in appendix C, and b2Q and b2i are given by 
equations (C14) in appendix C. 

Potential energy of the external forces .- The potential energy of the external forces 
applied to the ring is 



where F a , F q, F r , and are the stress and moment resultants in the shell which 

are functions of the deformations of the shell and where l\, I 2 , Z3, and Z4 contain the 
contribution of the applied external stress and moment resultants at the shell boundary. 

Potential energy of the pressurizing gas .- The equations developed up to this point 
in the analysis are applicable to rings of both open and closed cross section. It is also 
useful to extend the analysis to incorporate a pressurized ring of arbitrary closed cross 
section. For such a cross section, the change in the potential energy of the pressurizing 
gas must be accounted for in determining the stiffness of the ring. The change in the 
potential energy of the pressurizing gas has been derived in reference 5 by utilizing the 
following assumptions: 

(1) The pressurized ring is considered to be a membrane in that the local bending 
stiffness of the walls are neglected. 

(2) The internal pressure in the ring is constant during deformation. 

In the notation of the present paper and consistent with these assumptions, the potential 
energy of the pressurizing gas of the ring in terms of the shell variables and ring trans- 
verse shear strains is 


U D — pAj 



fc 


2r|L 


+ (A') + (B’r + 2CA' 


’ r c de 


(6) 


where A, B, and C are given by equations (3). 
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Shell Boundary Conditions 

Derivation of general boundary conditions .- The equilibrium equations of the ring 
and subsequently the boundary conditions of the shell can be obtained by variation of the 
total potential energy of the ring. From equations (4), (5), and (6), the result is 



+ EI z e 0 6^^ + EI xz e 0 6t]^ 


b 20 (p” + b 21 u '7( b 2O 6 0 M + b 21 6u ’’/ 


+ A r Gjy yz 6y yz + y yx 6y yx ) + ^ r s ^b 20 4>' + b 21 u ')( b 2O 6 0 r 



corresponding to arbitrary variations of the shell middle-surface displacements and rota- 
tion 6u, Sv, 6w, and 8<p and arbitrary variations of the ring shear strains 6y yx an d 
6y yZ . Equation (7) can be rearranged to isolate the generalized forces in the ring which 
correspond to the generalized displacements u, v, w, cp, y yx , and y yz . For equilib- 
rium of the ring, equation (7) vanishes for arbitrary virtual displacements, leading to the 
following equations which must be satisfied: 


rF a ± r c Q x = l 1 

(Or 

6u = 0) 

(8a) 

r¥ e ± r c Q 2 = l 2 

(Or 

6v = 0) 

(8b) 

rF r ± r c Q 3 = l 3 

(Or 

6w = 0) 

(8c) 

rM v ± r c Q 4 = Z 4 

(Or 

6<p = 0) 

(8d) 
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q 5 = o 

(° r 6y yx - 

(8e) 

0 

II 

CO 

<y 

(Or 6y yz = 0) 

(8f) 


where Q-^Qg,. . .,Qg are listed in appendix D. The choice of the plus or minus signs 
in equations (8) depends on the origin chosen for the shell coordinates. Equations (8a) 
to (8d) become the general nonlinear boundary conditions of an arbitrary shell of revo- 
lution with an attached ring provided equations (8e) and (8f) are used to eliminate the 
shear strains y yz and y^. Before carrying out this elimination, however, it is 
appropriate to express equations (8) in a linearized form. 

Linearization of the boundary conditions .- Many practical problems of aerospace 
shell structures involve small asymmetric deviations away from an axisymmetric state 
of equilibrium. The ring boundary conditions will be written herein in a form applicable 
to such problems wherein the asymmetric linear deviation may be due to either bifurca- 
tion buckling or the application of an asymmetric load. To this end, the stress and 
moment resultants ^F a , Fg, F r , and M^; displacements, rotation, and shear strains 
^u, v, w, <p, y yx , and y yz j ; and external loads ^ are expressed as the sum of two 
parts; for example, 


F a = F a,A + F a,B ( 9 ) 

where the term with subscript A describes an initial axisymmetric equilibrium state 
and the term with subscript B is an incremental change in the term away from this 
equilibrium state. 

The linearized set of boundary conditions for the incremental equilibrium state 
(subscript B) is obtained by substituting equations typified by equation (9) into equa- 
tions (8), subtracting out the initial state (subscript A), and neglecting terms nonlinear 
with respect to subscript B variables. The hoop force, moments, and strain in the ring 
for the initial and incremental states were also obtained and are listed in appendix D. The 
linearized subscript B variables are then taken to vary harmonically in the circumfer- 
ential direction in the following manner. The quantities Ug, Wg, <£g, F a? g, F r> B> 

B> ^1 B’ ^3 B’ ant * ^4,B vary in a cosine distribution; for example, 

U B = u b^ cos n0 (10) 

and v B , y yZjB , y^ g, F^g , and ig, B vary sinusoidally; for example, 

Vg = v^sin n6 (11) 
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The boundary-condition equations for the incremental state can then be put in matrix form as 


to 


) 


F w 

F a,B 

F (n) 

F 0,B 

F r"L 


0 

J 


G 11,A + G 11,B I G 12,A + G 12,B 


± r. 


g 22,A + g 22,B 


G 13,A + G 13,B 


g 23,A + g 23,B 
g 33,A + g 33,B 


Symmetric 


G 14,A + G 14,B 


g 24,A + g 24,B 
g 34,A + g 34,B 


G 44,A + G 44,B 


l 1 “ 

G 15,A + G 15,B g 16,A + g 16,B 

uW 

B 

1 1 __ 

1 g 25,A + g 25,B J g 26,A + G 26,B 

v w 

g 35,A + g 35,B ' g 36,A + G 36,B 
1 1 

w|.) 

G 45,A + G 45,B G 46,A + °46,B 

1 1 

" | “ - -- -- -- -J 

; g 55,A + g 55,B ', G 56,A + G 56,B 
1 1 

y(n) 

yz,B 

, i G 66,A + G 66,B 

y(n) 

yx,B 

v J 




r ^ 
,(n) 

,(n) 

*2,B 

M 

‘ 3,6 


,(n) 

4,B 


> ( 12 ) 


0 

v j 


where the G.. A and G.. „ are given in equations (El) and (E2), respectively, in appendix E. The subscripts A 

1J,A 1J,C 

and B refer to terms arising from the initial and incremental stress states, respectively. 


For special conditions, these general coefficients simplify. For example, for a ring which has no internal 
pressure, has no restraint of warping, has the shear center and centroid coincident, and is attached at the centroid 
to the middle surface of the shell, equations (12) reduce to 



i 


If it is further assumed that I xz = 0, equations (13) become, after some rearrangement, 



mi ii mi tin linn i 



Equations (14) are convenient for use in a subsequent section to compare results from the present theory with well- 
known results in the literature for in-plane and out-of-plane ring buckling. 

Final form of linearized boundary conditions .- To obtain the stiffness coefficients for the ring in the form of 

equations (1), the shear strains and _ must be eliminated from equations (12). The last two equations 

of equations (12) give 


•\ 


(n) _ ^56 \ 16 U B ^ + ^26 V B^ + ^36 W B^ + ^46^ B/ 

y yz,B ~ 


G IG U W + G + G + G d>( n ) 

o6\ 15 U B + *55 V B + Lr 35 W B + U 45^B , 


G 55 G 66 “ G 56 


> U5) 


(n) _ ^56^15 U B^ + ^25 V B^ + ^35 W B^ + ^45^B^) ~ ^55^16 U B^ + ^26 V B^ + ^36 W B^ + ^46 ^b! 


yx,B 


rv rw /w O 

g 55 g 66 " g 56 


where 


Substitution of equations (15) into the first four equations of equations (12) leads to the 
desired stiffness coefficients in the form of equations (1) where 




G 11 = r c G 11 + 


~ G 15 G 66 + 2G 15 G 16 G 56 " G 16 G 55| 


G 55 G 66 “ G 56 


g 12 = r c 


G 13 ~ r c 


G 14 = r c 


G l 5 (* 

G,n + ~^1 


G 56 G 26 _ G 66 G 25, 


'12 


G 13 + 


G 14 + 


i) + ^16^ 


G 56 G 25 " G 55 G 26/ 


~ ^ ~2 
g 55 g 66 " g 56 


G 15\ G 56 G 36 ” G 66 G 35 ) + G 16V G 56 G 35 " G 55 G 36, 


G 55 G 66 " G 56 


G 1s( G 56 G 46 ~ G 66 G 45) + G 16( ( 


G 56 G 45 " G 55 G 46/ 


G 22 ~ r c G 22 + 


G 55 G 66 " G l6 

~ G 25 G 66 + 2G 25 G 26 G 56 ' G 26 G 55j 

iv /v IV 2 

G 55 G 66 “ G 56 


G 23 ” r c 


G 24 = r c 


G 23 + 


G 24 + 


G 25( G 56 G 36 " G 66 G 3s) + G 2s( G 56 G 35 ~ G 55 G 36 


n > ( 17 ) 


G 55 G 66 " G !6 


G 2s( G 56 G 46 ~ G 66 G 4s) + G 26^ 


G 56 G 45 " G 55 G 46/ 


G 33 " r c G 33 + 


G 55 G 66 " G 56 

~ G 35 G 66 + 2G 35 G 36 G 56 ~ G 36 G 55) 
G 55 G 66 ' G 56 


G 34 ~ r c 


G 44 = r c 


J 34 


i ^35^ 


G 56 G 46 " G 66 G 45, 


i) + G 36^ 


G 56 G 45 " G 55 G 46 ; 


G 55 G 66 “ G 56 


/ ivO /v IV IV /V O I**' 

j P 6 r« . no r« r« n" p 


G 44 + 


~ G 45 G 66 + 2G 45 G 46 G 56 ~ G 46 G 55) 
G 55 G 66 ‘ G 56 


If transverse shear strains in the ring can be neglected, equations (17) still define the 
stiffness coefficients of equations (1) where now 
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^15 = ^16 = ^25 = ^26 = ^35 = ^36 " ^45 " ^46 _ 0 

The form of equations (1) is appropriate for computer programs if the shell behavior is 
defined in cylindrical coordinates. The transformation required for shells defined in the 
intrinsic shell coordinates is given in appendix F. 

DISCUSSION OF RESULTS 
Boundary-Condition Equations 

Equations (1) are the general boundary conditions for an arbitrary shell of revolu- 
tion with an elastic ring of general cross section rigidly attached at the shell boundary. 
These equations as developed in this paper can be used for the stress analysis of sym- 
metrically or asymmetrically loaded shells with a symmetrically prestressed ring 
attached at the shell boundary. The terms represent the external loads, if any, 
applied at the shell boundary. 

When the equations are applied to stability problems, the external loading on the 
shell must be symmetric so that the prebuckling stresses and deformations of both the 
shell and the ring are also symmetric. (This restriction, of course, does not imply that 
the buckling deformations must be symmetric.) The prestress terms in the stiffness 
matrix are now obtained from the axisymmetric prebuckling solution. Also, for buckling 
problems, the l± are always zero if the external loads on the ring are zero. 

An exception to the previous discussion must be noted for the special case of the 
pressurized ring. For the pressurized ring, the boundary conditions given by equa- 
tions (1) are valid only for axisymmetric stress analysis and symmetric or asymmetric 
buckling analysis. Equations (1) are not valid for asymmetric stress analysis of a pres- 
surized ring because of the way in which the internal pressure terms were handled in the 
perturbation procedure. (When nonhomogeneous terms appear on the right-hand sides of 
eqs. (8), eqs. (10) and (11) are not valid solutions for n > 0.) 

Buckling Solutions for a Ring 

To check the validity of the ring stiffness coefficients developed in this paper, 
selected ring buckling calculations were carried out for a ring of doubly symmetric cross 
section subjected to a constant radial line loading q. Equations (14) are in a convenient 
form to consider both in-plane and out-of-plane buckling. The axisymmetric prestress 
deformations are obtained from equations (14) by setting n = 0, letting = -r c q, and 

replacing the subscript B with A so that 
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i ii i nmi in 1 1 mw ■■ 


W A = — 


- 9 

q r c 


EI X 
EA r + — 

<*■ ft 


= 0 


y yz,A 

<p A = o 


y yx,A = 0 


j 


(18) 


The other two prestress deformations v A and vanish from axisymmetry and lack 

of rigid body motion, respectively. The prestress strain, hoop force, and moments may 
be determined by substitution of equations (18) and (3) into equations (D4) which yields 


q r r 


-\ 


'o,A 


EAr 


EE 


n a = 


M. 


x,A 


EI x q 


EL, 

EA r + — ^ 
r c 


M z,A = 0 


(19) 


Since in the formulation of the buckling problem the incremental loads are zero, the deter- 
minant of the coefficients of the displacements of equations (14) must equal zero. Because 
M A is equal to zero, in-plane buckling and out-of-plane buckling are uncoupled. 

Z,A 


In -plane buckling . - The in-plane buckling solutions are obtained by setting the 
determinant of the coefficients of the in-plane variables of equations (14) equal to zero. 

A first-order approximation to the critical load is obtained by retaining only the linear 
terms in q in the expansion for the determinant. For the lowest buckling mode, n = 2, 
this expansion yields 


nr^ 

qr c _ 4 _ 4 

EIx 1 + 2 T + 3 i + 24 1 

ST+1 S S T + l" 


( 20 ) 
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where 


T = 
S = 


EA r rf 

EI X l 
GA r r2 f 


2EI X J 


( 21 ) 


The critical load as determined by equation (20) agrees with the result given by 
Ratzersdorfer (ref. 6) when S — 


Out-of -plane buckling . - The out -of -plane buckling solutions are obtained by setting 
the determinant of the coefficients of the out-of-plane variables of equations (14) equal to 
zero. Again a first-order approximation to the critical load is obtained by retaining only 
the linear terms in q in the expansion. For the lowest buckling mode, n = 2, this 
expansion yields 


where 




ei 7 


4 + l + iL 


/4 + I + A 


\p fS T + IV ip fS/ T + 1 


2iip 


(— 4f + 6g) + 


2f 2 S| 





Ek I 

EI Z > 


EI X + EI Z 



J 


( 22 ) 


(23) 


When T — 00 and S — equation (22) reduces to the classical solution contained in 
reference 7. 


Application to the Stress Analysis of a Conical Shell 

To illustrate the application of the ring theory to a shell problem, the boundary- 
condition matrices were used to solve a representative shell stress -analysis problem. 
The problem chosen is a conical shell subjected to a normal pressure loading q which 
is constant along the meridian and varies harmonically around the circumference of the 
shell (fig. 2(a)). Thus, q = q n cos n0 where q n is a constant. This loading is appro- 
priate because a general asymmetric load on a shell of revolution can be expanded into a 
series of such components. 
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The shell is simply supported at the small end (that is, all displacements and the 
meridional bending-moment resultant vanish) and contains a Z-section ring rigidly 
attached at the large end. The properties of the ring are shown in figure 2(b). 



Detail A 

(a) Geometry of conical shell with detail of various ring attachment positions. 



(b) Properties of ring. 


A r = 0. 92 in 2 (5.94 cm 2 ) 

I u =0.446 in 4 (18.56 cm 4 ) 
l 22 = 0.551 in 4 (22.93 cm 4 ) 
l 12 * -0.392 in 4 (-16.32 cm 4 ) 
J- 0.012 in 4 (0.50 cm 4 ) 
r * 0. 148 in 6 (39. 74 cm 6 ) 

E - 30 x 10 6 psi (207 GN/m 2 ) 

G ■= 10 x 10 6 psi (68.9 GN/m 2 ) 


Figure 2.- Geometry of conical shell and edge-stiffening ring. 
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where r 0 is the radius of the middle surface of the conical shell at the large end. 
Poisson’s ratio was taken to be 0.3 in the calculations and the shell and ring were made 
of the same material. The amplitudes of meridional moment resultant, meridional stress 
resultant, and normal displacement as functions of distance from the stiffened edge (up to 
0.6 of s m ax) are shown in figures 3, 4, and 5 for the Fourier indices n = 3 and 8. The 
effects of eccentricity are shown in the figures for the ring centroid attached to the shell 
(z 0 = 0) and for a leg of the ring attached to the outside of the shell (zq = 0.925 inch 



0 .1 .2 .3 .4 .5 .6 0 .1 .2 .3 .4 .5 .6 

s,s max s,s max 

(a) n = 3. (b) n = 8. 

Figure 3.- Nondimensional meridional moment resultant m^ as a function of distance from stiffened edge of conical shell 

which is subjected to external pressure. 
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Results are also shown in the figures for the shell with the same boundary condi- 
tions at the small end and with boundary conditions denoted as "modified simple support" 
at the large end. The modified-simple-support boundary conditions are (see fig. 2(a)) as 
follows: Radial motion is restrained, movement in the direction of the axis of the cone is 
unrestrained, meridional moment vanishes, and circumferential displacement vanishes. 
These boundary conditions could be used to approximate the contribution of a stiff ring to 
shell behavior. 

Results were also obtained for the same ring- stiffened shell problems by using the 
ring boundary conditions contained in reference 1. The calculations based on the ring 
theory of reference 1 agree with the present calculations to within 5 percent and conse- 
quently are not shown in the figures. Further comparisons of results obtained with the 
two theories are given in reference 5. 

Plots for n = 3 in figures 3, 4, and 5 show that taking into account the eccentricity 
of the ring can cause appreciable changes in the behavior of the shell. Ring eccentricity, 
however, was not important for n = 8. The results also show that the modified-simple- 
support boundary conditions for the shell are poor approximations to the behavior of the 
ring-stiffened shell. The approximation for stress resultants appears to be better than 
that for displacements and moments. 

CONCLUDING REMARKS 

The contributions of ring stiffness to the boundary conditions for asymmetric 
behavior of a prestressed shell of revolution have been derived for an elastic ring of 
arbitrary cross section rigidly attached to the shell boundary. The ring behavior includes 
the effects of shear deformation, restraint of warping, ring torsion, out-of-plane bending, 
internal pressure, eccentricity, and axisymmetric prestress. The results for the ring 
stiffness matrix are presented in a form that is well suited for inclusion in a computer 
program for analysis of shells of revolution. The ring stiffness effects are applied to the 
sample case of a conical shell stiffened at one end by a ring and subjected to an asym- 
metric load. The results show that it is important to include properly the ring effect and 
ring eccentricity in order to determine accurately the stresses in the shell. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., December 9, 1969. 
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APPENDIX A 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 

The International System of Units (SI) was adopted by the Eleventh General Confer- 
ence on Weights and Measures, Paris, 1960 (ref. 3). Factors for converting the U.S. 
Customary Units used herein to the International System of Units are given in the following 
table: 


Physical quantity 

U.S. Customary 
Unit 

Conversion 

factor 

(*) 

SI Unit 
(**) 

Area 

in^ 

6.452 x 10- 4 

meters 2 (m2) 

Length 

in. 

2.54 X 10“ 2 

meters (m) 

Moment of inertia . . . 

in 4 

4.162 x 10" 7 

meters 4 (m 4 ) 

Stress 

psi = lbf/in^ 

6.895 x 10 3 

newtons/meter 2 (N/m2) 

Warping constant . . . 

in6 

2.685 x 10-10 

meters^ (m6) 


^Multiply value given in U.S. Customary Unit by conversion factor to obtain equiv- 
alent value in SI Unit. 

j|c j|( 

Prefixes to indicate multiples of SI Units are as follows: 


Prefix 

Multiple 

giga (G) 

109 

centi (c) 

10-2 
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APPENDIX B 


DERIVATION OF THE KINEMATIC RELATIONS FOR A CIRCULAR RING 
ATTACHED TO A SHELL OF REVOLUTION 

The kinematic relations are derived for a circular ring of arbitrary cross section 
undergoing shear deformations and rigidly attached to the surface of a sheU of revolution. 
These relations are then expressed in terms of the middle-surface shell displacements 
and rotation in a cylindrical coordinate system by requiring compatibility of the ring and 
shell at the point of attachment. (A more detailed analysis may be found in ref. 5.) 


Displacements of a Point in the Cross Section of a Circular Ring 

The coordinate system to be used for the ring (fig. Bl) is as follows: a coordinate 
x normal to the plane of the ring, a coordinate y along the shear-center axis of the 
ring so that y = r s 0, and a coordinate z measured radially inward from the shear cen- 
ter. The subscript s refers to the shear center. The displacement vector at any point 
in the ring is made up of the displacement 
components u a , v a , and w a in the direc- 
tions x, y, and z, respectively. The 
rotation /3 is assumed to be positive clock- 
wise about the y-axis through the shear cen- 
ters. (The positive rotation vector is tangent 
to the ring circumference in the plane of the 
ring and in the direction of increasing 6.) 

In defining /3, the thickness -compression 
strains e x and e z are assumed to be zero so that the projection of the deformed ring 
cross section into the xz -plane remains invariant during deformation. Note that this 
assumption does not preclude a warping displacement normal to the xz-plane. Consistent 
with this assumption, the shear strain e xz is taken to be zero throughout this analysis. 

The displacements u a , v a , and w a for moderately small rotations are, to first 
order, 

u a (z , 6) = u s ( 6) + z/3( Q) (B la) 



Figure B1.- Ring coordinate system. 


v a (x,z,0) = 


r s - z 
r s 


v s (0) - W^(0) - u' s (6) - zy yz (0) - xy yx (0) 


r s 


w a (x,£>) = w s (6>) - x/3(0) 


(Bib) 

(Blc) 
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APPENDIX B — Continued 

where the last two terms for v a in equation (Bib) are the effects of transverse shear 
deformations of the ring. (See ref. 5.) Note that equations (Bl) do not include a contri- 
bution from warping of the cross section. This contribution is considered separately in 
appendix C. 

For a thin shell, the displacements in a cylindrical coordinate system (fig. B2) are 


u 



Figure B2.- Cylindrical coordinate system for arbitrary 
shell of revolution. 


Up(0,T],z) = u(0,7j) + a 5 Z0(0,?}) 


Vp(0,T7,z) = v(0,T7) + z^w^i?) 

+ a2U*(0,?7) + a 3 v(0,7jj] 

w p (0,77,z) = w(0,tj) + a 4 z4>(0,77) 

J 

where 


a l 


a 2 


-\ 


_ J_ 
R 0 

a 3 = " a l 


a 4 = ra 2 

dr /dri 

a 5 = " ra l 

r 


J 


(B2) 


(B3) 


and r is the radius of the shell middle surface and is the radius of curvature 
which generates the shell circumference. Here, up, Vp, and Wp are shell displace- 
ments in the axial, circumferential, and radial directions, respectively, at a distance z 
from the shell middle surface, and u, v, and w are the corresponding displacements 
at the middle surface. 


Compatibility of Ring and Shell Displacements at the Shell Boundary 


The ring displacements given by equations (Bl) are now expressed in terms of the 
shell variables given by equations (B2) by requiring the displacements and rotation of the 
ring and shell to be the same at the point of attachment. (See fig. B3.) The line of attach- 
ment of the ring to the shell is at 77 = Constant. Therefore, at the point of attachment, 
the following conditions must be satisfied: 

u a (0,z) = -u p (e,? 7 ,z) 


v a (0,x,z) = V p (0,7?,z)l 
w a (0,x) = -Wp(0,T7,z) 


(z = z sh> x = X, Z = z) (B4) 


/3(0) = 4>{e) J 
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APPENDIX B — Continued 



Here Zg^ is the distance measured normal from the shell middle surface to the point of 
attachment, and x and z are the axial and radial distances from the shear center of 
the ring to the point of attachment. By means of equations (B2) and (B4), the displace- 
ments at any point in the ring cross section u a , v a , and w a can be expressed in terms 
of the shell variables and ring shear strains. By grouping in powers of x - x c and 
z - z c , the result is as given by equations (2) and (3) in the body of the paper where 


bj = x - x c - a 4 z sh 
i>2 = z - Zq + &5 z sh 


b 3 * Tpni 1 + a 3 z sh) 



(B5a) 

(B5b) 

(B5c) 

(B5d) 

(B5e) 

(B5f) 

(B5g) 

(B5h) 

(B5i) 
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APPENDIX B - Concluded 


b 10 = 




b ll 




b 13 = 0 
b 14 = 0 


b 15 



- z - z 
z 


c 


b 16 = “ 



- 1 


b 17 


r 


c 


r s - z 


x 


- X 


c 


b 18 = -l 



(B5j) 

(B5k) 

(B51) 

(B5m) 

(B5n) 

(B5o) 

(B5p) 

(B5q) 

(B5r) 

(B5s) 


In equations (B5), x c and z c are the axial and radial distances, respectively, from the 
shear center to the centroid of the ring, and r c is the radius to the centroid (fig. B3). 
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APPENDIX C 


STRAIN ENERGY FOR A CIRCULAR RING IN TERMS OF SHELL VARIABLES 


Strain Energy of Bending and Extension 

The nonlinear strain-displacement relation for the circumferential strain at any 
point in the ring in terms of the ring displacements (ref. 9) is, in the notation of the pres- 
ent paper, 



By substituting equations (2) into equation (Cl), the circumferential strain of the 
ring is expressed in terms of the shell variables and ring shear strains. If higher than 
second-order terms in z and x are neglected, the result is 

9 o 

e y = e Q + (z - z c )f?i + (x - x c )7 ] 2 + (z - z c ) V 3 + ( x “ x c) V 4 + ( x “ x c)( z “ z c)v^ ( C2 ) 


where 



(C3a) 

(C3b) 

(C3c) 

(C3d) 

(C3e) 

(C3f) 
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APPENDIX C - Continued 


and where 


\|xl = C’ - A 

A = (B f ) 2 + (C + A') 2 
B = 2B , 0’ - 2 -5-(C + A’)/ 


C 


-2(C + A , )l-~- 



D = (*•)- + 


2 . D 2 


r i 





j 


(C4) 


In equations (C4) the variables A, B, C, D, and F are given by equations (3). 

The strain energy for bending and extension of the ring due to the circumferential 
strain e y is 


= I L e ? dv - f L < l dA >-[ r c - <* - zcHm 


(C5) 


where V is the volume of ring material. Equation (C5) is expressed in terms of the 
shell variables by using equation (C2). Integration over the cross section of the ring 
yields 



h[V 2 + 2e o T ?4/ + 2Ix Z e 0 T7 5 + 


e o^ 


r c de 


(G6) 
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where use has been made of the following standard definitions of moments of inertia and 
area: 




™ x^jdAr — 0 



- Z JdAr- = 0 


> 



dA r - I z 


(C7) 



Equation (C6) can be expressed in terms of the ring hoop force and moments. First, 


N = 



M x = 
M z - 



(C8) 


where N is the ring hoop force and M x and M z are the ring bending moments about 
the x- and z-axes, respectively. Then, equations (C8) are integrated by substitution of 
equation (C2) into equations (C8) and by making use of equations (C7). The result is 


N - E^A r e 0 + I x 7?3 + Iz 7 ?4 + Ixz V 5 J 
M x = E^IxTjj + I XZ 77 2 ) 

M z = E^I z 772 + ^xz 7 ?^ 


(C9) 
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Finally, equations (C9) can be incorporated into equation (C6) to give the strain energy 
of the ring resulting from extension and bending as 


2ir/2Ne f 


A e 2 - — — M 

r6 ° t:t x 


“5r 7? i +- sr 77 2 i 


The first two terms of the integral in equation (CIO) represent the energy associated with 
stretching of the ring centroid line, the last two terms represent the energy associated 
with the in-plane and out-of-plane bending, and the third term represents a coupling of 
the bending and extensional behavior caused by ring curvature. 


Strain Energy of Restraint of Warping 

A rigorous development for the strain energy of restraint of warping would require 
a knowledge of the distribution of the warping displacements over the cross section of the 
ring and would constitute an extremely detailed and complicated analysis. Since restraint 
of warping is not a primary effect in ring behavior, a first approximation to restraint of 
warping was felt to be adequate for defining shell boundary conditions. This first approx- 
imation of the strain energy of restraint of warping is in the same form as that for a 
straight beam (ref. 7); namely, 

U 2 = f 2 V) 2 r c d0 (Cll) 

2r s r c J 0 L 

where £ is the twist of the ring with respect to the shear center and T is the warping 
constant of straight-beam theory. Values of r for cross sections of various shapes can 
be found in numerous structural handbooks. It might be noted that reference 10 also 
obtained the expression for strain energy of restraint of warping of a ring by making cer- 
tain simplifying assumptions on ring behavior. 

The twist of the ring (ref. 7) is expressed in the notation of the present paper by 


u 'j (C12) 

which is expressed in terms of the shell variables by use of equations (2), (3), (B4), 
and (B5). The result is 



5 =b 20 <*>' +b 21 u’ 


(C13) 
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where 




b 20 »^(l-b 9 ) 



(C14) 


Substituting equation (C13) into equation (Cll) gives the strain energy of restraint of 
warping in terms of the shell variables; namely, 




2r s r c ^0 


(Cl 5) 


Strain Energy of Shear 

In this analysis, the strain energy of shear is 


u 3 = § ( e yz + e yx)£c - ( z - z c) dAr <*0 


(C16) 


where V is the volume of ring material. It is convenient to integrate equation (C16) by 
separating the total shear strains into two parts — one part due to bending which is a 
function only of 6 and another part due to torsion which is a function of 6, x, and z. 
The result is 


e yz — “Vyz(^) Vyz(^>‘^>z)l 


eyx = - r yx + ryx( e ’ x > z ) 


(C17) 


where the y yz and y^ are the shear terms resulting from bending (the same as those 
used in eqs. (Bl)) and the y yz and y yx are the remaining shear terms resulting from 
torsion. By substitution of equations (C17) into equation (C16), the strain energy of shear 
becomes 


U 3 = 


A r G 


n2n 

J 0 


/ + y 

yx r yz. 


r c d 0 


+ § C f-2y, r „y, 7 , - 2 y v 


2 J-r 


yz'yz 


'yx^yx + V + V dV ( C18 ) 


where V is the volume of ring material. The first integral on the right-hand side of 
equation (C 18) represents the strain energy due to bending whereas the second integral 
represents the strain energy due to torsion plus a coupling of torsional and bending shear. 
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Clearly, these coupling terms are zero if the torsional strains are zero. Hence, if the 
coupling effect is neglected, the last integral reduces to the strain energy of pure torsion. 
As in the case of restraint of warping, a rigorous development of this part of the strain 
energy would necessitate a detailed study in itself. It is felt that this detailed study would 
not add materially to the results of the analysis and, consequently, the following simpli- 
fied expression, similar to that used in straight-beam theory (ref. 7), is used for the 
strain energy of torsion: 

“torsion <« 9 ) 


so that the total strain energy of shear now becomes 


2 2 \ , GJ r s 

V + v) r ° de+ T^ 


ClW + b 21 u’) 2 r c dfl (C20) 


where J is the torsional constant of straight -beam theory. 

In summary, the total strain energy for a ring of arbitrary cross section rigidly 
attached to a shell of revolution is given by the sum of equations (CIO), (C15), and (C20). 
The result is given by equation (4) in the body of the paper. 
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APPENDIX D 


LIST OF FUNCTIONS 


The following is a partial list of functions which are referred to in the main body of the paper and are included 
for completeness: 


«l-3 


M, 






*c r s 


~ B ' ■ ~ D + Pi3^i" x x 


M y - 


EI Y e, 


•x c o 


+ — M" 


( b 12^ ” b 10 A )( Mz 


EL 


L xz 

r c e °y 


r c 


y r c r s r 2 r c/ X ° 


lf£b 2 l(V>"" + b 2l“"’) 


- f r s b 2l( b 2O0" + b 21 u ") - ^(B" - b 10 «’) ♦ ^ B fe (Ae 0 )’ + ± 


r2 


p 13 A + b 12 § " 2b 12 £) EI xz e o 


(Dla) 




M. 


EI x e 0 \ b0 EI xz eQ\ j /P17 P 4 


•c / r 


~2 A ( M z " 


+ -y— $ - — DEI xz e 0 
r 2 r s r s j xz 0 


P 4 


— bo 


2 AEI x* e o " pA r ~2 ^ 


(Dlb) 


r c r§ 


\ r C / 




r c\ r c/ r c 


(p 54 .!iD)k-5ie 


£1 

r c 


AM z -^e, 


/p 18 ^ 
v r c 


CO 

C 71 


“§ D ) EI x e o ! 




(Die) 
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CO 

CT> 


Q 4 -~N--L 

1 c r f 


(-b 2 B- + P 2 ?)(n - 




♦!%-*) + !i i, 
r c\ r c/ r c x r| 


P2 


El, 


[p 8 B'-b 2 <K + P 9 S-^DHM x -^e 0 l 


+ — Ml' + -^ M„ - -i- 
r„ z r„ z „2 


(?i« { ' 



EIxz 


M 7 — e„ 

z r Oi 


£20 4 .£9 D+ p 0' + ^lBj EI x e o 

V*c *s *c / 


+ ^EI z (Ae 0 


p 9 A + p i6|;r " T f 


\ r S r c 


EIxze 0 > + b 20 (b 20 ^’”’ + b 21 u"") - SI r s b 2o (b 2o 0” + b 21 u") 

ISC c 




(Did) 


«5 = ^ - ^(N - ^ ^ + M D + p 15 |) ( Mx - 3S>) - ^ Mi ♦ i(b 18 | - 

c 



EIxz 


e °) + { ' ^ D ) EIxeo ' 7T EIzAe ° ' ^( p is A + it D - r a bl8£ ) EI - e ° + ArG V 

C ' ' c c\ ' 


M, 


- P A r|-y yx +_ |r 5) 


(Die) 


Q 6 


-^M^ + ^-^D + p 14 ^M x --^j.-^A(M z - — ej + 


T l 


■c V 


EI x e o\ b 15 \ L EI XZ \ .1 /P 2 2 


r C' C 


€ 


- ^ d) EI x e 0 - ^ AEI xz e 0 + A r Gy yz - pA r /-y yz + ^J 
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where 


Pl = b 5 - 1 

p 2 = b l + b 7 

P3 = 0 

P4 = -b4 
n b5 “ 1 

P 5 “ 

^7 

P6 = b 8 + r^ 

P7=7" 

r c 

- 2b 2 

P8 = 1 - 7T 


P9 =b 8 


, »(»>i - w) 


p 10 


r c r s 


Pll = “ 


r s 

c 


n XC 

P12 = 


X C b 10 

Pl3-^ri + TT 
b i5 - r s 

p 14 = — 

b 17 - x c 
p 15 = r c 

P16 = b 9 - 1 
p 17 = -b 4 

bs - 1 
r c 

2x c bm 
c r s A c 


Pl8 = 


p 19 r 


> (D2) 


p 20 = 2b 8 


, 3 ( b l + b 7) 


9 3b 2 

P21 = 2 - rT 


bis - 2r s 

P22 T r 


_ b i7 - 2x c 
P23 !> 


and 


| a C + A’ 




(D3) 


The terms bi,b2,- . .,b 2 j are given by equations (B5) and (C14); the terms A, B, C, 
D, and F are given by equations (3). 
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The strain, hoop force, and moments in the ring for the initial equilibrium state are 


a A 

e °,A = - rT 


"\ 


M E 
n A = TZ 


-A A Ar + |U A ^ 


m x,A - ^(" A A ^ + ^A^zj 


m z,A - ^f A A T~ + *M«) 


J 


(D4) 


The hoop force and moments in the ring for the incremental equilibrium state are 


n B = K 1 v b + K 2 w B + K3U B + K 4^B + K 5 w B + K 6^B + K 7Vyz,B + K 8 r yx,B 


M 


x,B = K 9 v^ + Kjo w b + K n u” + K 12 ^ + K 13 w b + K 14 0 B + K 15 y^ B + K 16 ^ xB \ 


M 


z,B = k 17 v b + K 18 w b + k 19 u b + K 20^b + K 21 W B + K 22^B + K 23Vy ZjB + K 24 > / ^x jB J 


(D5) 


where the Kj are given by 
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APPENDIX E 


o 


COEFFICIENTS OF RING STIFFNESS MATRIX 

A listing is provided of the coefficients for the ring stiffness matrix which appears in equations (12). 
coefficients which are a function of the initial state A are 


g 11,A = n2 


n A/ u 2 m x,A /, u 2 ou x c\ ^ m z,A EI x e o,A/ x c }\ 2EI xz e o A „ EI z e o,A 

^ 10 + 1 ) + ^. blo + 2 bl 0 _j + ^ bl0t ^_ tl+ ^ Xc+ _. 


— nbo / M„ a\ 

G 12,A = -4 N A b 10 + M x,Abll 


G 13,A = n2 i N A ^10^j) + ^ 


A ‘ r2 

1 n 


bgbiQ + 


( b 5 - *)' 


r c r s 


x c x s 


c 


(bl + b7)b 10 + b2 


^.A 


r c r s 


( b l + b7) + bgb 10 


bo n 

4- _ - 1 

M z,A 

( b l + b 7) 

T ± 

r c 

4 

L r s 


(b 9 - l) b io 


E I X e 0 ,A 

Xc 

r 4 

c 

r s 


(bi + b 7 + r c b 8 ) + b 2 - r c + EIx ^ e °? A bi + b 7 + r c bg + x c (b 9 - l) + EI y°? A ^b 9 - l) 

J r c r s J r c r s 


G 15,A = n 


N A ^ + ^( bl0 r s - r c b llbl5 ) - M ZA + Xc + 


r 

*c c 


r^r r’ 

c s * c 
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The prestress quantities appearing in equations (El) are defined in equations (D4) and the constants b^ are defined 
in equations (B5) and (C14). 

The coefficients which are a function of the incremental state B are: 
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APPENDIX F 


FORMULATION OF BOUNDARY CONDITIONS IN INTRINSIC COORDINATE SYSTEM 


The boundary conditions given by equations (1) are expressed in a cylindrical 
coordinate system and the shell variables are also in that coordinate system. However, 
the boundary conditions are usually difficult to use in this form in connection with shell 
programs that utilize normal and tangential displacements as variables. A more con- 
venient form is to express these boundary conditions in the cylindrical coordinate system 
in terms of the shell variables in the intrinsic coordinate system. 

From figure FI, it is noted that the stress resultants and middle-surface displace- 
ments in the cylindrical coordinate system are related to those in the intrinsic coordinate 



(a) Cylindrical coordinate system. <b) Intrinsic coordinate system. 

Figure FI.- Displacements and stress resultants in cylindrical 
and intrinsic coordinate systems. 


system by the following equations: 


F a = -Q^ cos a + N^ sin a 
F r = N^ cos a + sin a 


u = u sin a - w cos a 
w = w sin a + u cos a 


(FI) 


The circumferential displacement and rotation, as well as the stress and moment resul- 
tants Fg and M^, are the same in both coordinate systems. 
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By perturbing equations (FI), by subtracting out the initial equilibrium state (subscript A), and in addition by 
defining 

(n) 


%,B = cos ne 
Nt],B = n ^B cos ne 
ug = 4 n) cos n 0 




W B = Wg^cos n0 


(F2) 


J 


the boundary conditions in the cylindrical coordinate system in terms of the shell variables in the intrinsic coordi- 
nate system can be written as 
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In some instances, it may be advantageous to express the boundary conditions in the intrinsic coordinate sys- 
tem and in terms of the shell variables in that intrinsic coordinate system. By eliminating the coupling terms in the 
coefficient of the force vector in equations (F3), the following new form for the boundary conditions is obtained: 
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